WEAK APPROXIMATION OF STOCHASTIC PARTIAL DIFFERENTIAL 
EQUATIONS: THE NON LINEAR CASE 
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Abstract. We study the error of the Euler scheme applied to a stochastic partial differential 
equation. We prove that as it is often the case, the weak order of convergence is twice the 
strong order. A key ingredient in our proof is Malliavin calculus which enables us to get rid 
of the irregular terms of the error. We apply our method to the case a semilinear stochastic 
QO ' heat equation driven by a space-time white noise. 



< 



£ 



April 9, 2008 



< 

<-| '. 1. Introduction 

-i— > . 

When one considers a numerical scheme for a stochastic equation, two types of errors can 
be considered. The strong error measures the pathwise approximation of the true solution by 
a numerical one. This problem has been extensively studied in finite dimension for stochastic 
differential equations (see for instance [2D], [2B], [57], [52]) and also more recently in infinite 
dimension for various types of stochastic partial differential equations (SPDEs) (see among 
others P, [4], [6j, [TU], [11], p], p], [H], p], [IB], p], [E], [22], [23], [29], [30], [34], 
[35J, |36j). Another way to measure the error is the so-called weak order of convergence of 
a numerical scheme which is concerned with the approximation of the law of the solution at 

■^J- ■ a fixed time. In many applications, this error is more relevant. Pioneering work by Milstein 

([24], [25]) and Talay ([33]) have been followed by many articles (see references in the books 
cited above). Very few works exist in the literature for the weak approximation of solution of 
SPDEs. A delayed stochastic differential equation has been studied in [3]. Weak order for a 
SPDE has been studied only recently in [7J, [8], [19]. In order to explain the novelty of the 

/\ • present article, let us focus on a specific example. 

We consider a stochastic nonlinear heat equation in a bounded interval I = (a, b) C R with 
Dirichlet boundary conditions and driven by a space-time white noise: 



(1.1) 



^ = *« + f{X) + a(X)r,, £ € I, t > 0, 
X(a,t) =X(b,t) = 0, t >0, 

{ x(t,o) = x(0,£ei. 
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Where / and a are smooth Lipschitz functions from R to R. 

We introduce the classical abstract framework extensively used in the book |5j. We set 
H = L 2 (I), A = %, D(A) = H 2 (I) n Hl(I), W is a cylindrical Wiener process so that the 
space-time white noise is mathematically represented as the time derivative of W. We set 
f(x)(£) = f(x(C)), x £ H and define a : H^C(H) by a(x)h(£) = a(x(£))h(£), x,h£H. We 
then rewrite (II. ip as 

( , „, J dX = {AX + f(X))dt + a(X)dW, 

[l - Z) \ X(0) = x. 

It is well known that this equation has a unique solution. We investigate the error committed 
when approximating this solution by the solution of the Euler scheme 

(1 3) f X k+1 -X k = At (AX k+1 + f(X k )) + a(X k ) (W((k + l)At) - W(kAt)) , 

\ X = x, 

where At = T/N, N e N, T > 0. 

The study of the weak error aims to prove bounds of the type: 

\E(cp(X(nAt))) - R(tp(X n )\ < cAt 5 , 

with a constant c which may depend on <p, x, N and on the various parameter in the equation. 
Also ip is assumed to be a smooth function on H. If such a bound is true, we say that the 
scheme has weak order 5. In comparison, the strong error is given by M(\(X(nAt)) — X n \) or 
E(sup n=0 N \(X(nAt)) — X n \). Clearly, if the scheme has strong order 5 then it has weak 

order 8 > 6. Indeed, the test functions (p are Lipschitz. In general, it is expected that the weak 
order is larger than the the strong order. 

In the case of the Euler scheme applied to a stochastic differential equation, it is well known 
that the strong order is 1/2 whereas the weak order is 1 (see [32 J. The classical proof of 
this uses the Kolmogorov equation associated to the stochastic equation. The main difficulty 
to generalize this proof to the infinite dimensional equation (|1.2p is that this Kolmogorov 
equation is then a partial differential equation with an infinite number of variables and involving 
unbounded operators (see (|3.6|) below). The delayed stochastic differential equation studied 
in [3] is an infinite dimensional problem but since the equation does not contain differential 
operators the Kolmogorov equation is simpler to study. In [19|, a SPDE similar to (|1.2p is 
considered but very particular test functions (p are used. They are allowed to depend only 
on finite dimensional projections of the unknown and the bound of the weak error involves 
a constant which strongly depends on the dimension. In [7], [8], the Kolmogorov equation is 
not used directly. A change of variable is used in order to simplify it. In [7], the stochastic 
nonlinear Schrodinger equation is considered and the fact that the linear Schrodinger equation 
generates an invertible group is used in an essential way. This is obviously wrong for the heat 
equation considered here. The same change of unknown works in the case of a linear equation 
with additive noise as shown in [8] but there it is used that the solution can be written down 
explicitly. We have not been able to generalize this idea to the non linear equation considered 
here. 

We use in fact the original method developed by Talay in the finite dimensional case. The 
weak error is decomposed thanks to the Kolmogorov equations on each time step. Each term 
represents the error between the solution of the Kolmogorov equation on one time step and the 
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approximation given by the numerical solution. Due to the presence of unbounded operators, 
this apparently requires a lot of smoothness on the numerical solution. The main idea here 
is to observe that the non smooth part of the solutions of (J1.2D and ()1.3f) are contained in a 
stochastic integral. We get rid of this stochastic integral thanks to Malliavin calculus and an 
integration by part. We are thus able to prove that as expected the weak order is twice the 
strong order without artificial assumption except from a technical one on a. We restrict our 
presentation to the abstract equation above, a nonlinear heat equation driven by a space-time 
white noise. However, our method is general and can be used for more general equations as 
will be shown in future articles. Also, we only consider a semi-discretization in time. A full 
discretization will treated in forthcoming works. 

Note that the method developed here does allow to recover the result of [Sj. Indeed, in the 
Euler scheme f)1.3[) . the linear term is fully implicit and we cannot consider a scheme where it 
is partially implicit such as the theta-scheme considered in [8] . Note also that the proof below 
are much more complicated than in [8] and [7j- 

Malliavin calculus has already been used for the numerical analysis of stochastic equations. 
In [2], it is used to prove an expansion of the error of the Euler scheme for a stochastic 
differential equation under minimal assumptions on the test functions ip. This is a completely 
different idea and the Malliavin calculus is used completely differently. It is not clear that 
such ideas could be used for a SPDE. In a different spirit, Malliavin calculus is used in |31j to 
analyse adaptive schemes for the weak approximation of stochastic differential equations. 

Our method is much closer to the method developped in [21]. There, the Malliavin calculus 
is also used to get rid of a stochastic integral which appears when writting down the weak 
error. However, it is done in a global way and the error is not decomposed as in the present 
article. A fundamental feature of Kohatsu-Higa's method is that the Kolmogorov equation is 
not used so that more general stochastic equation can be can considered. The solution does 
not need to be markovian. However, no SPDE have been considered with this method. 



2. Preliminaries and main result 

We consider the following stochastic partial differential equation written in an abstract form 
in a Hilbert space H with norm | • | and inner product (•, •): 

( dX = (AX + f(X))dt + a(X)dW, 
{ ' \X(0) = x, 

where the unknown X is a random process on a probability space (Q, J-, IP) depending on t > 
and on the initial data x E H. The operator A is a negative self-adjoint operator on H with 
domain D(A) and has a compact inverse. We assume that 

(2.2) Tr((-,4)- a ) < oo, for all a > 1/2. 

We define classically the domain D((—A)P), /3 € R, of fractional powers of A and set 

\x\p = \(-Afx\, xeD((-Af). 
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The nonlinear function / takes values in H and is assumed to be C 3 with bounded derivatives 
up to order 3. We denote by Lf a constant such that for x,y £ H 

\f(x)\<L f (\x\ + l), 
(2-3) \f(x)-f(y)\<L f \x-y\, 

\f'(x)-f'(y)\ c(H) <L f \x-y\. 

The noise is written in terms of a cylindrical Wiener process W on H (see [5]) associated 
to a filtration (J r t)t>o- The nonlinear mapping acting on the noise maps H onto C(H), it is 
also assumed to be C 3 with bounded derivatives up to order 3. We denote by L a a constant 
satisfying 

/ 24) W(x)\c(H) < MM + !)> 

W(x) - a(y)\c(H) < L a\x - y\. 
We need a stronger assumption on this mapping, we require 

(2.5) \a"(x) ■ (h, h)\ c{H)) < L a \h\ 2 _ 1/4 , x£H, h £ H. 

Note that this implies a strong restriction on a. (See Remark 12.31 below for some comments 
on this assumptions). 

Recall that the cylindrical Wiener process can be written as 

where here and in the following (e^gjj is any orthonormal basis of H and (Pi)^^ is an 
associated sequence of independent brownian motions. This series does not converge in H but 
in any larger Hilbert space U such that the embedding H C U is Hilbert-Schmidt. Similarly, 
given a linear operator $ from H to a possibly different Hilbert space K, the Wiener process 
$W = Yle&$ 0£® e £ i s wen defined in K provided $ € £2{H,K), the space of Hilbert-Schmidt 
operators from H to K. (See the definition just below). 

Recall also that the stochastic integral L ^(s)dW(s) is defined as an element of K provided 
that \& is an adapted process with values in C2{H,K) such that f Q \ip(s)\ 2 c / H K \ds < oo a.s. 
(see [5]). 

If L £ C(H) is a nuclear operator, Tr(L) denotes the trace of the operator L, i.e. 

Tr(L) = y^(Lej, ej) < +oo. 

i>l 

It is well known that the previous definition does not depend on the choice of the Hilbertian 
basis. Moreover, the following properties hold for L nuclear and M bounded 

(2.6) Tr(LM) = Tr(ML), 
and, if L is also positive, 

(2.7) Tr(LM)<Tr(L)\\M\\ c(H) . 

Hilbert-Schmidt operators play also an important role. An operator L £ £-{H) is Hilbert- 
Schmidt if L*L is a nuclear operator on H. We denote by £-2(H) the space of such operators. 
It is a Hilbert space for the norm 

||L|| £2(H) = (Tr(L*L)) 1 / 2 = (Tr(LL*)) 1 / 2 . 
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It is classical that if L G C 2 (H), M G C(H), N G £{H) then NLM G £ 2 (#) and 

(2-8) \\NLM\\c 2i H) < ||^|U( ff) ||L|| £2(jff) ||M|| Ajff) . 

See [5], appendix C, or [9j for more details on nuclear and Hilbert-Schmidt operators. Note 
that (E2D implies that (-A)'? is Hilbert-Schmidt for any f3 > 1/4. 

Our assumptions imply that for any x G H, there exists a unique solution X(t) to equation 
to (|2.ip (see for instance [5], chapter 7). In the sequel, we often recall the dependence of the 
solution on the initial data by using the notation X(t, x). 

We approximate equation (|2.ip by an implicit Euler schemes. Let At = ^ > be a time 
step, we define the sequence (X k ) k =o ... N by 



/ 2 g\ J X k+ i = S At X k + AtSAtf(Xk) + V AtS A t^(X k )Xk+l, 

1 -X"o = x. 

We have set Xfc+i = (PF((fc + l)At) - W{kAt))/\/~At. The operators 5 A t is defined by 

S At = (I-AtA)-\ 

This is the classical fully implicit Euler scheme. It will be convenient to use the integral form 
of (I2TI) 



(2.10) X(t) = S(t)x + f S(t- s)f(X(s))ds + [ S(t- s)a(X(s))dW(s), t > 0, 

Jo Jo 

where S(t) = e is the semigroup generated by A. Similarly, (12. 9p can be rewritten as 



fe-i 



fc-i 



(2.11) X k = S k At x + AtY,S k A -/f(X e ) + s/AtY,S A -M X l)xi+x- 

1=0 1=0 

It will be convenient in the following to use the notation: 

t k = kAt, k = 0,...,N. 

The following inequalities are classical and easily proved using the spectral decomposition of 
A 



(2.12) 
(2.13) 

(2.14) 



(-AfS k At < ct, 13 , k>l, /?G[0,1]. 

(-AfS(t) < ct' 13 , t > 0, p > 0. 

C(H) ~ 



-AfS M <cAr p , 13 e [0,1]. 

C(H) 



(2.15) 



-Ay 13 (I - S At ) <cAt?, /3G[0,1]. 

C(H) 



Note that in (|2.1ip and (|2.9p . the noise term makes sense in H. Indeed, by (|2.14p . ()2.2p and 
(12. 8p . we know that S A t is a Hilbert-Schmidt operator on H. 
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We are interested in the approximation of the law of the solution of (12. ip . More precisely, 
we wish to prove an estimate on the error committed when approximating M((p(X{T,x))) by 
M((p(Xpf(x))). The function <p is a smooth function on H. 

In all the article, we use the notation Dp(x) for the differential of a C 1 function on H at 
the point x. If ip : i? i— > K, where K is another Hilbert space, Dip(x) G £(-ff, K) the space 
of continuous linear operator from H to K. When K = R, we identify the differential with 
the gradient thanks to Riesz identification theorem. We use the same notation and have the 
identity for x,h <E H: 

D(p(x).h= (Dp(x),h). 

Similarly, if ip G C 2 (H, R), D 2 ip(x) is a bilinear operator from H x H to R and can be identified 
with a linear operator on H through the identity: 

D 2 p(x).(h, k) = (D 2 ip(x)h, k), x,h,k£ H. 

Sometimes, we also use the notations ip', p" instead of Dp or D 2 ip. 

Given two Banach spaces K\ and K2, we denote by \\-\\k the norm on C^(Ki, K2), the space 
of k times continuously differentiable mapping from K\ to K2 with derivatives bounded up to 
order k. 

We use Malliavin calculus in the course of the proof. We now recall the basic definitions. (See 
|28j). Given a smooth real valued function F on H n and rpi,...,ip n £ L 2 (0, T, H), the Malliavin 
derivative of the smooth random variable F(f Q (ipi(s),dW(s)), . . . , f Q (ip n (s), dW(s))) at time 
s in the direction h € H is given by 



D h a 



/-•( / (Ms),dW(s)),..., [ (Ms),dW(s)) 



twff 



{iPl{s),dW(8)),...,J (Ms),dW(s)) } !r ; (,)./;) 



We also define the process DF by (DF(s),h) = D^F. It can be shown that -D defines a 
closable operator with values in L 2 (Q,L 2 (0,T,H)) and we denote by D 1 ' 2 the closure of the 
set of smooth random variables as above for the topology defined by the norm 



|F|| D i, 2 = (e(\F\ 2 )+E([ \D s F\ 2 ds\ 



1/2 



We define similarly the Malliavin derivative of random variables taking values in H. If G = 
Y,ieN F i e i G L 2 (Q,H) where Fi G B 1 ' 2 for all i G N and J] ieN J'q \D s Fi\ 2 ds < 00, we set 
D h a G = T,teN D s F iei, D S G = Y,i&i D sFie t . We define B 1 ' 2 (iJ) as the set of such random 
variables. 

When h = e m , we write D £m = D m . 

The chain rule is valid and given u G C£(R), F G B 1,2 then n(F) G B 1 ' 2 and D(u(F)) = 
u'{F)DF. Also if G = J2im F i e i G D 1,2 (H) and u G C^^R) then u(G) G B 1 ' 2 and 
D(u(G)) = Du(G).DG = (Du,DG), or equivalently D*(u(G)) = 52 ieN d l uD%F i = (Du,D%G). 

Note that as already mentionned, we identify the differential of a function in C 1 (-ff, R) with 
its gradient. 
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For F G B 1 ' 2 and V G L 2 (VL x [0,T];JT) such that ij){t) G B 1 ' 2 for all t G [0,T] and 
J" J q |L> s f/>(t)| 2 <i,s(it < oo, we have the integration by part formula: 

¥.(f J (ij(s),dW(s))] =Elf (D s F,^(s))ds 

where the stochastic integral is a Skohorod integral which is in fact defined by duality. In this 
article, we only need to consider the Skohorod integral of adapted processes in which case it 
corresponds with the Ito integral. Moreover, the integration by part formula above holds for 
F G B 1 ' 2 and i\) G L 2 (S7 x [Q,T];H) when ip is an adapted process. Recall that if F is T t 
measurable then D S F = for s > t. 

We will often use the following form of the integration by part formula whose proof is left 
to the reader. 

Lemma 2.1. Let F G B 1 ' 2 (iJ), u G C%(H) and rp G L 2 (Q x [0,T],£ 2 (H)) be an adapted 
process then 

e(du(F)-J i/;(s)dW(s)\ = ^(^f D 2 u(F)-{D™F^{s)e m )ds\ 



E ( f Tr (^*(s)D 2 u(F)D s F) dsj . 



Also we remark that this Lemma remains valid if u is not assumed to be bounded but only 
u G C 2 (H) provided the expectations and the integral above are well defined. This is easily 
seen by approximation of u by bounded functions. 

We now state our main result. 

Theorem 2.2. Assume that f and a are C^ functions from H to H and C(H) and that a 
satisfies (|2.4j) . then for any x G H, T > 0, e > 0, the Euler Scheme (12. 9p satisfies the following 
weak error estimate 

\E&(X(T,x))) -%(I W ))| < C(T,\^\ c! ,\x\,e)At 1 / 2 - £ , ^^(ff). 

Remark 2.3. Assumption (|2.4p is quite restrictive. It is void for an additive noise or a noise 
of the form BX dW where B is a linear operator from H to C(H). Otherwise, it implies that 
the noise is a perturbation of such noise. An example of a noise satisfying this is 

a{x) = Bx + d{{-A)- 1 ' 4 x) 

where B G C{H) and a : H^£(H) is a C 3 function with derivatives bounded up to order 3. 
This assumption is crucial in our proof. It is used in essential way in Lemma \4-5\ which is 
used at many points of the proof. 

Apart from this point, our result is optimal. If the noise is assumed to satisfied some non 
degeneracy assumptions, the smothness assumption on the test function ip can be weakened. 
This will be investigated in a future work. 

In all the article, C or c denote constants which may depend on A, f, a, Q or T but not 
on At. Their value may change from one line to another. The initial data x is fixed and the 
constant may also depend on |x|. Note also that we assume that At < 1, we could also assume 
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At < Aio for some Ato > 0. In this case, the different constants would depend on Ato- Finally, 
e is a small positive number. 

3. Proof of the main result 

The proof uses different tools from stochastic calculus such as ltd formula, Kolmogorov 
equations, Malliavin calculus. Sometimes, it may be very lengthy and technical to justify rig- 
orously their use in infinite dimension. We avoid these tedious justifications by using Gakerkin 
approximations. We replace equation (|2.ip by the finite dimensional stochastic equation 

dX m = (AX m + f m (X m ))dt + a m (X m )dW, X m (0) = P m 

where P m is the eigenprojector on the m first eignevectors of A, f m (x) = P m f(x), a m (x) = 
P m a(x)P m . It is not difficult to prove that X m converges to X in various senses. 

Similarly, we replace the discrete unknown X k by a finite dimensional sequence defined in 
an obvious way. 

We prove the result for these finite dimensional objects with constants that do not depend 
on the dimension m. It is then easy to deduce the result for our infinite dimensional equation. 

In order to lighten the notation, we omit to explicit the dependence on m below and write 
X, /, a instead of X m , f m , a m . 

Step 1: We first define a continuous interpolation of the discrete unknown. 
We rewrite (12.91) as follows: 



X k+1 =X k + A At X k + S At f(X k )ds + / S At a{X k )dW(s) 

Jt k Jt k 

where A A t = S A tA. Note that A A t is in fact a Yosida regularization of A and is a bounded 
operator: 

(3.1) \A At \ C {H) < cAr 1 . 
It is then natural to define X on [0, T] by 

(3.2) X{t) = X k + f A At X k + S At f(X k )ds + I S At a{X k )dW(s), t G [t k ,t k+1 ). 

Jt k Jt k 

Clearly, X is a continuous and adapted process. Given a smooth function G on [0, T] x H, ltd 
formula implies for t € [t k ,t k+ i) (see [5]): 

ft JIQ 

G(t,X(t)) =G(t k ,X(t k ))+ —{ S ,X{s)) + L k>At G(s,X{s))ds 

(3.3) 

+ f (DG(s,X(s)),a(X k )dW(s)). 
Jt k 

Where for ^ G C 2 (H,R) 

L k , At ip(x) = ^Tt {(S At a(X k )) (S At a(X k ))* D 2 i>(x)} + (A At X k + S At f{X k ),D^(x)). 

Step 2: Decomposition of the error. 
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Let us define 
(3.4) u(t, x) = E(ip(X(t, x))), t € [0, T\. 

Then the weak error at time T is equal to 

u(T, x) - E(<p(X N )) = E(u(T, x) - u(0, X N ) 

(3 ' 5) = £E(«(T-t fc ,* fc )-u(r-tfc + i,* fc+ i)). 

fc=0 

It is well known that u is a solution to the forward Kolmogorov equation: 

du 

— (t,x) =Lu(t,x) 

(3.6) 

= -Tt{a(x)a*(x)D 2 u(t, x)} + (Ax + f(x),Du(t, x)). 

Therefore, ltd formula (|3.3p implies 
E(u(T - t k+1 ,X k+1 )) = E(u(T - t k , X k )) + E / L kAt u(T - t, X(t)) - Lu(T - t, X(t))dt. 

The first term in (|3.5p will be treated separately and we decompose the error as follows 

N-X 
(3.7) u(T, x) - E{^{X N )) = u{T, x) - E{u{T - At, Xj) + ^ a k + b k + c k . 

fe=l 
Where 

a k = E I L+1 (AX(t) - A At X k , Du(T - t, X(t)j\ dt, 

J~kk 

b k = E (f(X(t)) - S At f(X k ),Du(T - t, X(t))\ dt, 

Jtk 

c k = -E / Tr | [a(X(t))a*(X(t)) - (S At a(X k )) (S At a(X k ))*\ D 2 u{T - t, X(t))j dt. 

^ tk 

In the next steps, we estimate separately the different terms in (|3.7p . 

Step 3: Estimate of u(T, x) - E(u{T - At, Xi)). 
By the Markov property 

u{T, x) = E(<p(X(T, x))) = E(u(T - At, X(At))). 
Therefore, by Lemma 14.41 f° r an y e > 0, 

\u{T, x) - E{u(T - At, X x ))\ < c(T - Aty^^MUE (\X(At) - X x \_ l/2+s ) . 
Moreover 

X(At) - X x = (S(At) - S At )x + J Ai S(t - s)f(X(s, x))ds - AtS At f(x) 



+ Jo S(t - s)a(X(s, x))dW{s) - V r AiS At a(x)xi- 
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It is easy to prove that 



| ( _ A) -l/2 +£ {s{At) _ S M )\ C(H) < cAt 1 ^- 



Since (S(t))t>o is a contraction semigroup and |(— ^4) l ' 2+e • | < c| • |, we have by (|2.3p and 
Lemma 14.21 

At 

< AtL f E( sup \X(s,x)\ + l) <cAt(\x\ + l). 
i/2+e se[o,At] 



E 
Similarly 



S{t-s)f(X(s,x))ds 
o 



\AtS At f(x)\_ 1/2+£ <cAt(\x\ + l). 

We then have 

»At 



E\\J S(t - s)a(X(s, x))dW{s)\\ /2+£ 

= E ^K-AyV^Sit - s)a(X(s,x))\l 2{H) ds 

< E (J * |(-^)- 1/2+£ |^ (H) |5(t - s)|^ (if) |cr(X( S , x))|^ (fl) (fa 



and by ([22]), ([23]), Lemma [O] 

EM/ 5(t - s)a(X(s, x))dW{s)W /2+ A < cAt{\x\ + 1). 

Similarly 

E (\^/AtS A t(r(x)xi\ 2 ) < cAt(\x\ + 1). 
Gathering these estimate and using Cauchy-Schwartz inequality , we obtain 
(3.8) \u(T, x) - E(u(T - At, X x ))\ < c{T - At)~ 1/2+E Ar l/2+£ < cAt l ' 2 - £ 

where, as mentionned above, the constant is allowed to depend on T, x, (p, /, a . . . 

Step 4: Estimate of aj-, k > 1. 
We split a k as follows: 



1 i 2 

a k — o-k + a k 



with 

'■tfc+i 



a\ = E I"' 1 ({A- A At )X k ,Du(T - t, X(t))) dt, 






al = EJ hl (A(X(t) - X k ), Du{T - t, X(t)) 



dt. 

Note that A At - A = 6AtS At A 2 . By Lemma l4~4l below, we know that Du(T - t,X(t)) is in 
D((-Ay) for 7 < 1/2 and it is easy to see that X k belongs to D((-A) s ) for 5 < 1/4. It is 
impossible to compensate the presence of A 2 by such arguments. The idea is to recall (|2.1ip 
and to observe that the irregularity of X k is contained in the stochastic integral. We thus 
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further decompose al in three terms according to (|2.1ip . The first two terms are easy to treat. 
The third one involves the stochastic integral and is estimated thanks to Malliavin calculus. 
We set 

a]; 1 = -6AtE J k+1 (s At A 2 Si t x,Du(T-t,X(t))) dt, 



a 1 ,' 2 = -9 AfE 



*fe+i 



fe-i 



t'k 



SMA 2 AtY J Stt i f{Xt),Du{T-t 1 X{t)) dt, 



e=o 



1,3 



-9AtE 



r t k+1 J fe-i \ 

/ ' [s At A 2 VAiJ2S k ^ e cr( X t)Xi+i,Du(T-t,X(t)))dt, 



so that 



1 1,1 , 1.2 . 1,3 

a k = a k +a k + a k . 



By ([iniD , (f2"7L2|) and LemmaO we have for k = 1, . . . , N - 2 and e > 

(3.9) 



\a]: l \ <cAm 



f fc+1 ISAti-A^^lc^-A^Si.lc^K-A^^DuiT-^Xm Wt 



< cAt 1 /2-2e t -l +£ f tk+ \ T _ t) -{l/2-e) dL 

Jtu 



1,2 



The estimate of a k is similar. We have by (|2.3|) . (j2.12[) 



fc-l 



At(-A)^Y,S k At e f(X £ 



e=o 



fc-1 



'At 



<L f AtJ2\(-A) X - e Sl 

1=0 

fe-i 



C(H) 



(1**1 + 1) 



fcO 



Since 



fc-l 



IrpE 



*tY,& ^ £ ' lT 



e=o 



we deduce thanks to Lemma 14.41 and Lemma 14.11 



(3.10) 



\a k \ < cAt 



t+i 



S'At(-A) 



l/2+2e 



C(H) 



(T - i)- 1 / 2 ^^ 



< cAi 1 / 2 ^ // ft+1 (T - t)-W 2 ~ £ Ut 
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To treat a k ' , we first rewrite it in terms of a stochastic integral and then use Lemma 12.11 



4' 3 = 6AtE [ k+ V I " '■' ' 2 ^'-' 
Jt k \Jo 

tk + l ftk 



S At A z S K At ts a(X es )dW{s),Du(T - t,X(t)) dt 



= 6AtE / Tr l(T*(X is )S At A 2 S A ^ is D 2 u(T-t,X(t))D s X(t)\dsdt 

where £ s = [s/At] is the integer part of s/At. By the chain rule and (|3.2p . we have for 
s E [0,t k ], /i G H,t G [t fc ,ifc+i), 

iftyM r-jft -v' i / /i n^ v i c t'fv \ r^h v a„ i / c I *Uv \ r^h 



D h s X{t) = D h s X k + / A At D h s X k + S A tf'(X k ) ■ D h s X k ds + / S At {<j'(X k ) ■ D h s X k \ dW(s). 
Jt k Jt k 

Fro < 1/4, we have, by (I2TT4D . (12^) . ([2^1) . (1231) 



E 



/ S/^t (cr'{X k ) ■ D h s X k \ dW(s) 
Jti, 



E 



< E 



t k 



t k 



(~AfS At (a'{X k ) ■ D h s X k 



C 2 (H) 



ds 



(-A) 



-1/4-e 



C 2 (H) 



{ _ A) P+l/4+e SAt a i {Xk) . D h Xk 



C(H) 



C(H) 



ds 



<cAt l / 2 -^- 2£ E{\D h s X k \ 2 ). 

We then use (|3.ip . f|2.3|) to bound the other terms above and obtain thanks to Poincare in- 
equality 



(3.11) E(\D*X(t)\j) <cE[\D h s X k \^, 8€[0,t k ], te[t k ,t k+1 ). 

By Lemma 14.31 we obtain for (3 < 1/4 



E 



(-AfD s X(t) 



C(H) 



< ct 



-2/3 
fc-4" 



1.3 



We are now ready to conclude the estimate of a k ' . We choose e > and write thanks to ([2? 
(EUD, (El2|), Lemma H3] and (H2D 



|aj' 3 | < 9AtE 



*fc+l ftk 



t h JO 



W*(Xi 



\C{H) 



S At A l l 2+2£ {-Af^^S^ 

C(H) m 



C(H) 



{-A) l / 2 ~ £ l 2 D 2 u{T-t,X{t)){-A) l ' 2 - £ l 2 Trli-Ay 1 / 2 -^ 2 } (-AfD s X(t) 

< cAtE [ k+1 f At~ l / 2 ~ 2e r k l ^ £/2 {T - t)~ 1+ %l ea ds dt. 
Jt k Jo 



C(H) 



ds dt 
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-t k + -l+e/2. 



Since £ k t~_J s £/2 ds < § T £ / 2 , we deduce 

(3.12) |4' 3 | < cAt 1 / 2 " 26 / k+1 (T - t)~ 1+e dt. 

Jt k 

Gathering flUD, (l3TT0|) and (l3TT2|) . we obtain for k = 1, . . . , N - 1 
(3-13) 



|4| < cAt 1 ^ {t k l+£ + 1) (f k+ \T - t)- l+£ dt + l) . 



We now estimate a 2 . Let us set 



4 1 =E [ k+ \t- t k ) ( AA At X k , Du(T - t, X(t)) ) dt, 



tk 

tk 
tk+1 ft 



af = E / (t-t k )( AS At f(X k ), Du(T - t, X(t)) ) dt, 



a 2 / 3 =E I I [AS At (r(X k )dW(s),Du(T-t,X(t)))dt, 

Jt k Jtk v J 

so that thanks to (|3.2p . we have a\ = a k ' + a fc ' + a fe ' . The first term a k ' is similar to a\ 
above and is majorized in the same way 

(3.14) loj 1 ! < cAt l ' 2 ~ 2e (t- 1+e + 1) (f k+ \T-t)- 1+£ dt+l\ . 

for k = 1, . . . ,N — 1. The second one is not difficult to treat, we have using similar arguments 
as above 



(3.15) 



< cAtK-A^SMicmmim) / cr-r (1/2 ~ £) ^ 

KcAt 1 / 2 - 6 f tk+ \T-t)- {1/2 - £) dt 
Jtk 



for fe = 1, . . . 3 iV — 1. The estimate of aj requires the use of Lemma 12 .11 It implies 



a,' = E 



/ Tr {a*(X k )S At AD 2 u(T - t,X(t))D s X(t)\ ds dt. 



'tk 
Since, X k is Tt k measurable, we have from (|3.2|) 

(3.16) D s X{t) = S At a{X k ) s E (t k ,t k+1 ], t k <s<t<t k+1 . 



14 A. DEBUSSCHE 

It follows, thanks to flZ3D, (jZZHj) , ([Z2]) and Lemma H31 
a} 3 = E / (t - t fc )TV [a*(X fc )S , AtAD 2 «(T - i, X(£))S Ai a(X fe )} (it 



<cAtE/ kC^i^^ifirAtC-^^^^k^K-Aj^-^^^^.^^^^.^i^-e/ai 

xTr((-A)- 1 /2- £ /2 ) | ( _ Ar5At | £W | a(Xfc) | i:(H)dt 

< cAi l/2-3e/2 f tk+ \ T _ t )-l+e dt 
Jt k 
for fc = 1 , . . . , iV — 1 . Finally, we obtain 

(3.17) |4| < cM l ' 2 - 2£ {t k l+£ + l)(f k+1 (T - ty 1+£ dt + l\ 
for A; = 1, . . . , N — 1. Together with ()3.13|) this yields the estimate of a k 

\a k \ < cM l ' 2 - 2e {tl l+£ + 1) (f k+ \T - ty 1+£ dt + l) . 
It follows easily 

N-l 

(3.18) Yj l«fel < cAt 1 / 2 " 22 . 

fc=i 

Step 5: Estimate of b k . 

This term seems easier to treat since we do not have the unbounded operator A. However, 
since it involves the nonlinear term, we need to use Ito formula (|3.3p to control f{X(t)) — f(X k ), 
this introduces many terms. For some of them we again use Malliavin integration by parts. 

First, we get rid of S A t ■ We have thanks to (|2.15p . (|2.3p . Lemma PP1 and Lemma HTTJ 

b\ = E I L+1 ((I - S At ) f(X k ), Du(T - t, X(t)j\ dt 



< cE / (1 + \X k \)\(-A)-^ z+£ (I - S At )\c iH )\(-A)^ z - £ D U (T - t,X(t))\dt 

tk+l 



< cAt 1 / 2 -^ f k+1 (T - ty 1 / 2+£ dt 

Jtl. 



>t h 

for k = 0, . . . , N — 1. We now estimate 



b 2 k =bk-b\ 

= E I*' 1 (f(X(t) - f(X k ),Du(T - t,X(t)j\ dt 



'tk 



E / " ' 5D(/i(^(*) - fi(X k ))d iU (T - t, X(t))dt, 

Jtk „-,-lM 



tk i&N 
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where fi = (/, e{) and di = (D-, ej). We choose (e^).^ as the orthonormal basis of eigenvectors 
of A. By (pOI) . we have for i E N 



/i(X(t) = fi(X k ) + / ^Tr {(S At a(X k ))(S At a(X k )yD 2 f i (X(s))} ds 

f (A At X k + SAtfiX^DfiiXis))^) ds + j (Dfi(jt(s)), a{X k )) dW{i 



+ f 

tfc ' Jt k 

With obvious notations, this defines the decomposition 

,2 ,2,1 . ,2,2 . ,2,3 . ,2,4 

To treat the first term, we rewrite it as followsj: 



b 



2,1 



1 ftk+1 ft 

-E 



^Tr [(S At a(X k ))(S At a(X k ))*D 2 fi(X(s))} d lU {T - t,X(t))dsdt 



t k J tk jgf^ 

= I E /"* ' / Tr {(S At a(X k ))(S At a(X k )yA(s,t)}dsdt 
Jtk Jtk 

where A(s, t) £ C(H) is defined by 

(A(s, t)h, k) = Y, D 2 fi(X{s)).{h, k)d iU {T - t, X(t)) 



D 2 f(X(s)).(h, k),Du(T - t,X(t)) , h,keH. 



Obviously 



\As,t)\c( H )< D 2 f{X(s)) 



C 2 {HxH,H) 



Du(T-t,X(t)) 



where £ 2 (H x H,H) denotes the space of bilinear operators from H x H to H. By (|2.3p and 
Lemma 14.41 we deduce: 

\A(s,t)\ c{H) < c. 
Then, we write thanks to flZgp , (|2J4"|) . (pP]) . 

|Tc{(5 A to-(X fc ))(5Ata(Xfc))*A(a > t)}| 

< Tr ((-^l)- 1/2 - £ ) K-^+^a^) k(^)l^) |^(«,t)| Afl) 

^cAt-^-^l + lXfcl) 2 . 



We deduce by Lemma 14.11 
(3.19) 



bf < cAt 3 / 2 - £ . 



Recall that we in fact work with Galerkin approximations so that all sums below are finite sums. 
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The second term b k ' involves the same difficulty as a\ above. We rewrite it using p. lip . This 
gives 



b l' 2 = E 



+E 



ftk + l ft ( k—1 \ 

/ / Y[AAtSi t x + A At AtJ^S A - l f(X e ),Df i (X( S )))d i u(T-t,X(t))dsdt 

Jtk Jt k iGN y i=0 J 



«GN 



, t J k S^aiXi^IWirlDfAX'^))) <)MT - I- X{l))d*dl 



where, as above, £ T = [r/At]. The first term is bounded as follows, using Lemma 14.41 dZ2]), 
(12T2JL d2Hj), 

E f" /l^ [AAtSitX + AAtAt^Sfc'fiXelDfiiXis))) d iU (T - 1, X(t))dsdt 



'th J t k jgN 



1=0 



E 



/ fe+ Y ( D /(^( s ))- U At S k At x + A At AtJ2s k A ~ e f(X e )) ,Du(T-t,X(t))) dsdt 



e=o 
fe-i 



< cE r +1 f\{-AfS At \ c{H) ( \{-A) l ^S k At x\ + j; k-A) 1 "^ 



— e ok— I. 
At 



C{H) 



\f{X t )\ )dsdt 



<cAt 2 - £ (t k 1+e + l) 



(.2,2 



The second term of 6,' requires an integration by parts, we obtain 



E 



[^ I Y, ( AAt I h S^aiXt^mir^Df^Xis))] d % u{T - t,X(t))dsdt 
Jt k Jt k , cM V Jo J 



rtk+l rt 



E l I Y, ( A ^ / S A -^a(X eT )e m ,e j )d(5 m (r)d j f l (X(s))d i u(T-t,X(t))d S dt 

'th Jt k ■ ■ rM \ JO 



E 



ifcj-i rt rt 



t h JO 



53 ( j4 At^/ T ^(^Je m ,e i 



i,j,m,n£N 



d j>n fi(X(s)) [D™X(s),e n )d iU (T-t,X(t)) 



+d j f l (X(s))d i>n u(T - t,X(t)) D?X(t), e. 



drdsdt 



E I l f /* 53 Z} 2 /^(s))(A At< S^<TpQ T )e m , D?X(sj) 8 lU (T - t,X(t)) 



tk Jtk JO 



i,m£N 



jk-tr 



+ Bi(s, t)A At S A - tT a(Xe T )e m , D?X(t) drdsdt 



E 



tk+i ft rtk 



tk Jtk Jo 



J2^{(DrX(s)yD 2 f i (X(s))A At S A - t ^a(X iT )}d iU (T-t,X(t)) 

+Tr { (p r X{t)f B t ( S ,t)A At S k A / T a{X er )} drdsdt 
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where, for i G N, Bi(s,t) is defined by 

(Bi(a, t)g, h) = {Df t {X{s)),g) £ ^, n u(T - t, X(t))(/i, e n ), 5 , h € JT. 

neN 

The first term above is estimate as b k ' . For the second term, we write 
£ ieN (Bi(s, t)g, h) = D 2 u(T - t,X(t)) ■ (Df(X(s)) ■ g, h) 

= {D 2 u{T - t, X(t))h, Df{X{s)) ■ g), g,heH. 
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Therefore 






< 



C(H) 



Df(X(s)) D 2 u(T-t,X(t)) 



C(H) 



We deduce by Lemma SU (ETTT1) . <|2^|> . <|2131t . (1231) . (1X121) . (|23jh Lemma O and similar 
arguments as above 

E / " +1 / E (^* / "" S k A -^a(X u )dW(T),Dfi(X(s))] d iU (T-t,X(t))dsdt 
Jt k Jt k icn \ Jo J 



i&i 

< cAt^ 2 - £ . 
Therefore 

It is also easy to see that 



bf < cAt 3 ' 2 ^. 



b l' 3 = ]E 



I" l I y / (SAtf(X k ),Df i (X(s)))d i u(T-t,X(t))dsdt 

Jtu Jtk ,^u V J 



= E f" 1 f Du(T - t,X(t)) • (Df(X(s)) ■ S At f(X k )) dsdt 
Jt k Jt k v ' 

< cAt 2 . 

2 4 

It remains to estimate b k ' . We again integrate by parts the stochastic integral and obtain by 
Lemma 12.11 



bl A =E 



E 



~k + l 



t k Jt k j gN 

tk + 1 ft 



tk Jtk 

t k + l ft 



E 

'tfe, Jt k 

< cAt 3 / 2 ~ £ , 



Y^ (Dfi(X(a)), cj{X k )dW(s)) 8 lU (T - t,X(t))dt 

ieN 

Tr { (l> s X(t)) * D 2 u(T - t, X(t))Df(X(s))a(X k )} ds dt 
Tr i(T*(X k )S At D 2 u{T - t, X(t))Df(X(s))a(X k )\ ds dt 



thanks to (ETI51) . (123J) and ([2TT3D . 



is 
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We conclude this step by gathering the previous estimates. This enables us to write 

JV-l 

J2 \h\ < cte 1 ' 2 -' 



fe=i 



Step 6: Estimate of Cfe. 

Using the symmetry of Du, we introduce the decomposition of c k : 



c k = ~® 



1 ftk+i 






Tr { [a(X(t))a*(X(t)) - (S At a(X k )) (S At a(X k ))*] D 2 u(T - t, X(t))] 



dt 



+ 



+ 



^EJ kl Tt[(I- S At )a(X(t)) ((I - S At )a(X(t))y D 2 u(T - t,X(t))] 

Tr [s At a(X(t)) ((/ - S At )a(X(t))^j * D 2 u(T - t, X(t))} dt 
^J" 1 Tr {s At (a(X(t)) - a(X k )) (s At a{X{t)))* D 2 u{T - t, l(i))} 
^E f * _1 Tr [s At a(X k ) (s M a(X{t)) - a{X k )) * D 2 u(T - t, X(t))} dt 



dt 



dt 



cl + 4 + 4 + 4- 



The first two terms are easy to treat, we use similar arguments as in the previous steps and 
write thanks to (|2.7p . Lemma H~5| Lemma I3~T1 f|2.15j) 



I < cE f k "' TtU-A)- 1 ' 2+£ (I - S At )a(X(t)a*(X(t))(I - S Ai )(-^r 1/2+£ } (T - t) 



-l+2e 



dt 



<<^j kl Ty[(-A)- 1 ' 2+£ {I - S At )(I - 5 Ai )(-^)- 1/2+£ } (T - t)- 1 

ptk+i 



-l+2e 



dt 



< cAt 1 / 2 - 35 / (T - t)- l+£ dt 
Jt k 

The second term is similar, we have 



c 2 <E 



•fe+i 



A)-y^(i-s At 



C(H) 



^^o\ c{H) \(-Ars At \ c{H) Tr{(-A)-y^} 

dt 



{-Af/ 2 - £ D 2 u{T - t, X{t)){-A) l l 



2-e 



C(H) 



<cAt 1/2 ~ 3s f k ^(T-t)- 1+£ dt 

Jti. 
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The estimate of the next term is much more complicated. It is based on similar arguments as 
before but the computations are much longer. 
We use (I3.3P and obtain for h,k £ H: 

V 

Ith 

ft 



a(X(t)) - a{X k )) h,k) =\J Tr {{S At a(X k )) (S At a(X k ))* D 2 (a(-)h, k) (X(s))} dt 



+\ f (A At X k + S At f(X k ), D (a(-)h, k) (X(s))) dt 

= (Ah,k) + {Bh,k) + (Ch,k). 



Thus we may write 

4 =-E Tv[S At A[S At a(X(t))) D 2 u{T - t,X(t))j 

+^E f " ^ Tr [s At B (s At a(X(t))Y D 2 u(T - t, X(t))} « 

+^E I " ^ Tr {s At C (s At <j(X{t))) * D 2 u(T - t,X(t))} i 



3,1 . 3,2 . 3,3 
c k + c fc + c fc • 



Note that 



(Ah,k) = \l Y, (^"(X(s)).(S At a(X k )ee,S At a(X k )e e )) h,k) ds. 

By (J23I), for u,v£ H, 

( (cr"(X(s)).(u,v)j h,kj< L (T \u\_ 1 / 4 \v\_ 1 / 4 \h\ \k\ < c\u\ \v\ \h\ \k\. 

We deduce, thanks to ([22]), (|2T4|) . 

(Ah, k) < c/\t l ' 2 - £ (l + iXfel) 2 ^! |fc|, 

and 

\A\ c{H) <c/xt 1 l 2 ~ F ~(l + \X k \) 2 . 

Then, by Lemma KT\ Lemma li~5| (J2.14J) and again ([212]) 



sj 1 < cAt 1 / 2 ^ f k+ \T-t)- l / 2+£ dt. 
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The term c, ' involves the same difficulty as at and b k ' . We use (|2.1ip to replace X% by a sum 



of three terms: 

(Bh, k) = \ I (A At S k At x + AtA At V S k A /f(X e ) + f * A At S A -^a(X iT )dW(r) 
1 Jt h \ , ,-, -Jo 

D{a(-)h,k)(X{s)))ds 



fc-i 



e=o 



+ k—1 

\ j ( \a'(X(s)) ■ {A At S k At x + AtA At £ S k A /f(X 



£=0 

t h a . c k-e T 



+ J, k A At S k A ^a(X eT )dW(r)) 



h,k)ds 



= {B 1 h,k) + (B 2 h,k) + (B 3 h,k). 
We then write thanks to flZ3D, flZIZD and (ETUI) 

< c 



a'{X{s))-A At S k At x h,k\ds 
o'{X{s)) ■ A At S k At x \h\ \k\ds 



C{H) 



< cAt\A At S k At x\ \h\ \k\ 

< cAt l - £ t\r e \h\\k\. 



Similarly 



(B 2 h,k) < cAt l ~ £ \h\\k\. 
It follows, thanks to Lemma |45j (J2l4l) and ([22]) 



h J " ~* Tr {5 A t(Bi + B 2 ) (s At a(X{t)j) * D 2 u(T - t, X(t))} 



dt 



<cAt l - 6£ (tr £ + l) / (T-t 



•fe+i 



'dt. 



3,2 



The estimate of the part of c,' involving ,63 is very technical. As before, we get rid of the 
stochastic integral thanks to an integration by parts. This results in a supplementary trace 
term. In order to work with the double trace, we write everything in terms of the components 
of the operators and vectors. Given an operator G on H, we set G 1 ' 3 = (Gei,ej). We thus 
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write 

E i ^ Tr [s At B z (s At a{X{t))) * D 2 u{T - t, X(t))} dt 

Jtk 

ftk+l — . . . _ / _ \ rn.i 

= E J2 &3 3 <7 m ' 3 (X(t))){SAtD 2 u(T-t,X(t))S At ) dt 

tk i,j,meN 



Y, E r l f f " d r a^(X{s)) [A At S k ^a(Xe T )e n , e r ) d(3 n (r)a m ^ (X (t))) 

_ \m,i 

S At D 2 u(T-t,X(t))S At ) dsdt. 



It is important to recall here that in fact we work with finite dimensional approximations of 
the solutions so that all the above sums are finite. We now use the Malliavin integration by 
parts and obtain 

E i ^ Tr (s At B 3 (s At a(X(t))) * D 2 u{T - t, X(t))\ dt 
Jtk 

£ E f k+1 f f tk ^d r ^(X( 8 )) (D n T X{s),e p ) (A At S k A -^a(X iT )e n ,e r ) a m ^X(t))) 



i,j,m,n,r€N h tfc u pSN 



s M i^u{T-t,x{t))SAt 



+ Y J d r a^{X{s)) (A At S A -^a{X ir )e n ,e r ) d p a m ^(X(t))) (D?X{t),e p ) [S At D 2 u{T -t,X(t))S At 

pen 



ik-e T 



+d r a^(X(s)) [A At S^a(X, T )e n ,e r )a m ^(X{t))) [S At [D 6 u(T - t,X(t)) ■ D^X(s)) S At ) drdsdt 

= 1 + 11 + III. 
We then write 



tk+l ft ftk 



I= E E / / / D 2 ^ 3 \X(s)).^X(s),A At S k At ^a(X eT )e n )a m ^X(t)) 

S At D 2 u(T-t,X(t))S At ) drdsdt 



^E / / / D 2 u(T-t,X(t))- U 1 (s,r,k)e j ,S At a(X(t))e j ) drdsdt 



i'gn * Jtk jo 



ifc+1 ft ft 



E / ' L f I "" Tr{ a* (X(t))S At D 2 u(T -t,X{t))fa{8,T, k)\ drdsdt 
h k Jt k Jo l > 



22 



A. DEBUSSCHE 



where we have set 

M*,r,k)hi = ^S A t (D 2 o{X{ S )) ■ ^X(s),A At S k At ^a(X iT )e n ^ hi, hi G H. 

neN 

Let us define £«,/n,/i a by 

(Ps,huh2 u > v ) = [ SAt { d2(7 (X(s)) ■ (u,v)j hi,h 2 ) , u,v £ H. 

Then by (123]) 

l^a./ii.haUfH) < c\hi\ \h 2 \. 
We deduce by ([23]), (J52J, (12TT41 . ([2A2]) . ([3TT]) and LemmaSJ 

(<f>i(s,T,k)hi,h 2 ) =Tc{**(X tr )S ] £ eT AteE athlM D T X(s)} 

<\a*(Xac(H)TT(-Ar 1/2 -%-A)y 2+s S k A -^A At \ ciH) \^^^^ 

< cM- l / 2 - 2 % l+£ \hi\\h 2 \(l + \X k \) 

and by Lemma 14.14 Lemma 14.51 and (|2.2p 

"*fe+i 



/ < C At 1 /2-2 £tfc l+ £ /" fc+1 (T _ t) -V2-e dL 
Jtk 



Similarly, we may write 



nen 



•fe+i 



t rt 



//= E E / / / Tr\[[Da(X(s))-[A At S k A -^a(X eT )e 

(-Da(X(f)) • (.D™X(*))) *S A t-D 2 u(r-t,X(t))5AtUrds(it 
- dr ;: -E / / / \<t> 2 (T,s,t,k)\ L{H) (T-t)- x / 2+e dTdsdt 



U JO 



with 



Mr,s,t,k) = Y^ [(Da(X(s)) ■ (A At S A -^a(Xe T )e n ))} [(Da(X(t)) ■ (p?X(t) 

neN 

We use similar arguments to estimate its norm. For u,v £ H, we write 
(<f> 2 (T,s,t,k)u,v) 

= J2neN(\(Dv(X(t)) ■ (D?X(tj))Y u, \(Da(X(s)) ■ (A At S k A ^ a(X, T )e 



Tr {a* (X iT )S k At ^ A At a* v b u } 



with 



a v h 



S At (Da(X(s))-h 
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Kh 



Since 

we deduce 

and 

Finally 



S M [Da(X(t)) • (D h T X(t)^ 

\ a v\c(H) <c\v\, \K\c(H) < c\u\, 
\Mr,s,t,k)\c(H) < cTr{S k A -^A At } < cAr 1 ^^ . 

II < cAt 1 ' 2 ^ f k+ \T - t)- l / 2 - £ dt. 
Jt k 

III=\ I" l f [ " V Tr{~f n S At a(X(t))}dT ds dt 

* Jtk Jtk JO _^M 



where for u,v G H 



n<EK 



Y,DMT-t,X(t)) (d?X(8),u,S m (Da(X(s)) ■ (A At S k At ^a(X eT )e n 



n€N 



= TrU(u,v)(-A)- 1 / 2 - £ (-A) 2£ D?X(s)\ , 
and for hi,Ji2 £ H 
(K(u,v)h!M) = r> 3 u(r-t,X(t))-((-A) 1 /2-^ 1 ,i t ,5 A t (dct(X(s)) ■ [A At S k A -^a(X eT )h 2 

By Lemma 14.61 

Ku,v)\ C (H) < ciT-ty^+X^Ar^ul \v\. 
Therefore, by (jgSjI . ([3TT|) and Lemma E3 

£>««,«) < C ( T " i)~ 1/2+£ Ct £A ^ 3£ l n l M- 

nGN 

It follows 

hnlcw^ciT-tr^+H^At 

and by ([22]), (gUJ) 



-3e 



/// < cAt 1/2 - 4e i k+ \T - t)- l ' 2+£ dt. 
Jtk 



We can now conclude 

I 3,2 



< cAt l ' 2 - A£ {t k l+£ + 1)( / fc+1 (T - ty l/2+£ dt + 1). 

Finally, it is easy to check 

\C\ C{H) < cAt(l + \X k \) 
and 

I 3,3i 



< cAt l - 2£ ["^(T- t)- l ' 2+£ dt. 
Jtk 



24 A. DEBUSSCHE 

We deduce 



; ftk + 1 

< cAi 1 / 2 - 4 ^ 1+£ + 1) / ((T - ty l / 2 - £ + l)dt, 

Jt k 

d in exactly the same way, 

< cAt 1 / 2 - 46 ^ 1 ^ + 1) f^\(T-t)- 1+£ + l)dt. 

Jt k 

N-l 



and, since c 4 , is majorized in exactly the same way 

rtk+l 

*tk 

It follows 

N-l 



fc=l 

Step 7: Conclusion. 

It is now easy to gather all previous estimates in (|3.7p and deduce 

\u(T,x)-E(ip(X N ))\ <cAt l / 2 ~ i£ . 

Recall that all the above computations have been done on the Galerkin approximations of X 
and X k ). The constant c above does not depend on m so that we can easily let m— >oo in this 
estimate and obtain the result. 

4. Auxiliary Lemmas 

In this section, we state and prove technical Lemmas used in the preceeding section. Again, 
the various estimates used here could be difficult to justify rigorously on the infinite dimensional 
equation and we in fact work with Galerkin approximations. Taking the limit m^oo at the 
end of the proofs gives the results rigorously. 

The first two Lemmas are very classical and we state them without proof. 

Lemma 4.1. For any I £ N ; there exists a constant ci such that 

max E(\X k \ l )<ci(\x\ l + l). 

k=0,...,N 

Lemma 4.2. For any I € N, there exists a constant c~i such that 

supE(|X(t,x)|) <ci(\x\ l + 1). 

t 

Lemma 4.3. For any (3 E [0, 1/4), there exists a constant c such for k = 1, . . . ,N , s £ [0,tk], 
we have 

tf_ e E (\D*X k \f) < c\h\ 2 , heH. 

Proof: By (|2.1ip and the chain rule, we obtain the following formula for the Malliavin 
derivative of X k : 

fe-i 
D h s X k =S k A - e *a(X is )h + At Y, S k A - e f(X £ )-D<:x e 

e=e s +i 

fc-i 
+VA1 £ S k At e (a'(X e ) ■ D f s l X £ ) X e+i 
e=e s +i 
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for s G [0, tk] and h G H. 

By (|212]> . (j23]l . (j2lS]l . (J231), we deduce for e > 



fc-i 



e=e s +i 



k-t 

At 



fe-1 



-At ^ (-A)^VTO-^ 



£(H) |/'(^)|£ W 



£fl(«) 



I&X* 



< c(t-li\hf + L%[At £ t 



/3 
fc-£ 






£>?X* 



2 



It is now easy to use a discrete Gronwall Lemma and prove 



,2d 



iX+tJ^" SXA ) - r|/ ' 1 



□ 



Lemma 4.4. Lei y? G C^(il, R). For am/ /3 < 1/2, i/tere exists a constant c@ such that for 
t>0, x G H 

\Du(t,x)\/3 < cpt~~P\\(p\\i, 

where u is defined in (|3.4p . 



Proof: Differentiating (|3.4p . we obtain for h £ H: 

Du(t, x)-h = E (D<p(X(t, x)) • v h ' x (t)) 

where n h,x (t) is the solution of 

dr] h,x = ( Ar] h,x + f(x(t, x)) ■ r] h ' x ) dt + a'(X(t, x)) ■ r] h ' x dW, 

v h > x (o) = h. 

We rewrite this equation in the integral form 

r, h ' x (t) = S(t)h+ f S(t-s)f'(X{s,x))-r l h ' x {s)ds+ [ S(t-s)a'(X{s,x))-r ! h ' x {s)dW(s), t > 0. 
Jo Jo 

By ([13]), (pOj) . (I2TT2]) . we have for j/, & G H and a > 1/2: 
\S(t)a'(y).k\ <L a {-A)^ 2 



C 2 (H) 



(-A) a / 2 S(t) \k\ < ct~ a/2 \k\ 

C(H) ~ 
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Using (|2.3|) and then Cauchy-Schwarz inequality, we obtain 

h,x 



E 



V 



(*)f) <cr 2 e\h\l + LpU£\r, h >*(8)\<b) )+Ej\t-sr a \ri h >*(sf 



ds 



< c r 2f5 \h\ 2 _ p + c f E{\r] h ' x (s)\ 2 ds + E / (t - s)~ a 7] h ' x (s 
Jo Jo 



ds. 



sup t 2/3 E ( 
te[o,T] V 



v h ' x (t) <Ni/j. 



It is classical that this implies 
(4.1) 

We deduce 

\Du(t,x) -h\ < c\\(p\\it~P\h\-0. 
Taking the supremum over h yields the result. 



□ 



Lemma 4.5. Let tp G C 2 (H, M). For any /3, 7 < 1/2, i/iere exists a constant cp~ such that for 
t>0, x £ H 

\(-AfD 2 u(t,x)(-Ay\ c{H) < C/3i7 H^)|M| 2 , 

where u is defined in (|3.4j) . 



Proof: We use the same notations as in the proof of Lemma 14.41 We differentiate a second 
time (13.41) and obtain for h,k £ H: 



(4.2) D 2 u(t,x) ■ (h,k)=E(D 2 ip(X(t,x)) • fa** (*),»/*•*(*)) + D<p(X(t,x)) ■ C hAx (tj) 

where C, h ' k,x {t) is the solution of 

' d( h > k ' x = (AC, h > k ' x + f"(X(t, x)) ■ {rj h ' x (t),ri k ' x (t)) + f'(X(t, x)) ■ ( h > k > x (t)) dt 

+ (a"(X(t, x)) • {ri h > x (t),ri k > x (t)) + o-'(X(t, x)) ■ C, h > k ' x (t)) dW, 

^ ^M(o) =0. 
We rewrite this equation in the integral form 

C fc '*>*(i) = y"5(* - a) (/"(*(*, a:)) • (r?^( S ),^( S )) + /'(*(*, x)) • C^*)) tfa 

+ |* S(t - s) (a"(X(s, x)) ■ (v h ' x (s),v k ' x (s)) + a'(X(s, x)) ■ ( h > k > x (s)) dW(s), t > 0. 
Using similar argument as above and (12. 5D . we prove 



(4.3) 



rt \ 2 

E(|C fc ' fc,x (t)| 2 ) <cE( I \r) h ' x (s)\\r) k ' x (s)\ + \C h ' k ' x (s)\ds 



+cE At - s)- a (\v h > x (s)t 1/A \ri k > x ( S )£ l/4 + \C hAx (s)\ 2 ) ds, t > 0. 
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Proceeding as in Lemma \4A\ we have thanks to Burkholder inequality and then to Minkowsky 
inequality 



E 



(\v h ' x (t)\ 4 ) <cr^\h\i fi + cE((r\ v h ^(s)\ds^ j+EN r 

< c t-W\h\i + c (J* (E(|^( S )| 4 )) V4 ds 
+c(j\t-sr a {E(\r, h ' x (s 



(t-s)~ a rj h ' x (s) 



ds 



Lx/„\|4' 



Taking the square root of this inequality and using a generalized Gronwall Lemma, we deduce 



sup t 4/3 E 
te[o,T] 



V n ' x (t) )<c\h\ 



(4.4) 

Similarly, we have 

E(I^WIii/ 4 ) <ct l - ip \h\t + c¥.(j \ V h > x (s)\ds) +tE( / (/- ,)-"|,/'- J '!.s)!^/., 



4 / ,t \ 2 

ai„ft,a:/'.,\|2. 



1/4 



ds 



+c(f(t- s)~ a E (|T7 h ' x («) | 4 ) V2 ds] 

Therefore, by (jOl) . 

(4.5) E(|^(t)| 4 _ 1/4 ) < rf 1 "^/^ 

Plugging these inequalities and similar ones for r\ ,x in (|4.3p yields 



sup E 
te[o,T] 



C ft '*(t) J < clfcl^l^. 
The result follows easily using (|4.4p and this inequality in (|4.2p . 



D 



The following Lemma is proved thanks to similar arguments. 



Lemma 4.6. Let if 6 C^(H,M). For any (3 < 1/2, there exists a constant cq such that for 
t>0,xeH,hie D{{-Af), h 2 £H, h 3 £ H 

D 3 u{t,x) • {{-AfhiMM) < c/^IMIsM \h 2 \ M, 
where u is defined in (|3.4p . 
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